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Abstract. - It is shown that the magnetic energy of a quantum current, contrary to the 
classical case, is essentially negative. Since this result allows to escape a famous theorem by 
Bloch, it can be expected that, under appropriate conditions, the ground state of a quantum 
conductor may be characterized by a spontaneous current. A similar idea was suggested, many 
years ago, by Heisenberg and by Born and Cheng as a possible mechanism for superconductivity. 



Introduction. - It is well-known, from the study of the infra-red divergence problem in 
relativistic QED, that a charged particle cannot be separated from its classical field 
It is known also that the motion of a classical particle cannot be described correctly, if the 
interaction with the self-generated field is not taken into account properly py. 

In this paper we analyse the effects induced by the self-generated classical {i.e. coherent) 
field on a stationary current flowing in a macroscopic quantum conductor. In particular we 
will estimate the contribution to the total energy of the system due to the classical magnetic 
field and to its interaction with the current. The analysis will concern explicitly two geomet- 
rically different models. In both cases we will obtain the remarkable result that the magnetic 
energy associated to the current is essentially negative. A completely different result would 
be obtained for a current due to a system of classical charged particles. Such a difference de- 
ll" 1 ■ scends from the fact that, in classical mechanics, the kinetic energy of a particle is expressed, 
in a natural way, in terms of the velocity. We remark that the velocity is essentially a clas- 
sical concept. In quantum mechanics the velocity of an electron corresponds to the following 
q | operator 



X 



v ^ m^ 1 [-ihV + eA c (t, r)], (1) 



where A c is the classical vector potential. Owing to the presence of A c , which depends on 
the other particles of the system as well as on the external environment, the velocity v, given 
by eq.(l), is essentially a collective operator. On the other hand, the canonical momentum 
is connected directly to the particle wave-length, which in turn is involved by a boundary 
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condition like periodicity. For these reasons it looks more proper, in our opinion, to analyse 
a quantum many-body system in terms of the canonical momenta and not in terms of the 
velocities of the particles. 

As previously mentioned, in this paper we will obtain the result that the magnetic self- 
interaction energy of a quantum current is essentially negative. This result may have un- 
expected consequences concerning e.g. the ground state of a system of electrons subject to 
periodic boundary conditions. In fact, as it will be discussed in the sequel, it is conceiv- 
able that, under appropriate conditions, a physical situation corresponding to a spontaneous 
current may be favoured with respect to a current-less situation. We recall that attempts 
to explain superconductivity in terms of spontaneous currents have been made long ago by 
Heisenberg 0] and by Born and Cheng 5 . However such theories have been ruled out by a 
theorem by Bloch 6 , owing to the lack, in our opinion, of a proper analysis of the magnetic 
self-interaction. 

The Schrodinger equation for a quantum system interacting with its classical field is given 

by 

ihd t \i>) =H[A C ] (2) 

where 

H[A C ] = H M - JdrA c - (j p + i j A ). (3) 

We have assumed the Coulomb gauge. In eq.(3) Hm represents the Hamiltonian for the 
matter, including the Coulomb interaction among charged particles. In the interaction term 
the current-density operator j has been splitted into two contributions. The first one, which 
will be called canonical current, consists of 

j p (r) = -ihJ2 a {q a /2m a )[5{r - r a )V a + V a 5(r - r a )]. (4) 
The second contribution, involving the potential A c , consists of 

3A{r) = -J2 faa/ m «) n ^ r ) Mr), (5) 

* — J a 

where n a (r) = S(r — r a ) represents the density operator for the a-th particle. A third contri- 
bution due to the spin magnetization-density has been neglected, for the sake of simplicity. 
The classical potential A c is solution to the Maxwell equation 

V x B c - c- 2 t E ±c = e^c- 2 {i;\{j p +j A )\i;)^ (6) 
with the transverse electric field given by 

E ±c = -d t A c . (7) 

Both eq.(6) and eq.(7) can be obtained, in the framework of QED, by use of the coherent-state 
formalism [Jj. 

In this paper we limit ourselves to consider the interaction of the system with the classical 
field, neglecting any residual interaction with the quantum field (uncoherent photon emission 
and reabsorption) . In this approximation, a simple expression for the conserved total energy 
of the whole system (matter plus classical field) can be given, in the form: 

£= (i>\H[A c ]\rP) + (e /2) f air (\E ±C \ 2 + c 2 \B c \ 2 ). (8) 



F. Miglietta : Magnetic energy of a quantum current 



3 



It can be verified easily that £ is conserved. In fact, by differentiating both sides of eq.(8) 
with respect to time, we obtain, in the absence of emission of radiation, 



dt 



r 8 a 

- eo c*V x B C - e „e»V ■ (E, x flj, . -^j ds ■ (Jfc, x flj . 0, 

where eq.s (2) and (6) have been used. 

The second term in the r.h.s. of eq.(8) represents the (positive) energy of the classical 
field. By a simple calculation, this term can be cast in the form 

| J dr(\E ±c \ 2 + c 2 \B c \ 2 ) = ^ J drA c -{i>\{j p +j A )\i>)+e J dr\E ±c \ 2 (10) 
J dr \A C \ 2 + ie c 2 <j> ds ■ (A c x B c ), 



4 dt 2 



where eq.s (6) and (7) have been used. In the absence of a significant emission of coherent 
radiation the last term can be neglected, for a finite system. By use of this result we obtain 
from eq.(8) 

£ = {H M )-\ J drA c -(j p )+e Q J dr\E ±c \ 2 -^^ J dr\A c \ 2 . (11) 

We notice the absence of any contribution due to ja ■ We notice also the factor of 1/2, as 
well as the minus sign, in the second term in the r.h.s. of eq.(ll). In stationary conditions, 
eq.(ll) reads 

£ = £ M +£f = (H m ) -^JdrA c - (j p ), (12) 

where the Coulomb gauge is understood. The analysis of the physical content of eq.(12), in 
particular of the minus sign in the r.h.s., can bc simplified if the potential A c is expressed in 
terms of the canonical current j p and not in terms of the total current appearing in the r.h.s. 
of the Maxwell equation (6). This result will be obtained explicitly for two models. 

First model. - The first model consists of a hollow cylinder with a current flowing around 
it. Let L be the height, R the inner radius, a the thickness and let us assume, in order to 
simplify the calculations, a -C R <C L. Let us assume an uniform canonical current, flowing 
around the cylinder, given by 

(j p ) = —en(hp/m) — —en(h[i/mR), (13) 

where n is the (constant) electron density. We observe that fj, is related to the total orbital 
angular momentum UM of the electrons by the relation /i = MN" 1 , where = 2nRaLn 
represents the total number of electrons involved in the process. 

From eq.(6) we obtain the following expression for the magnetic flux 

$ B ~ -fte _1 /i[l + ( 2 /7 2 a-R)] _1 - -he' 1 ^ (14) 

where 7 is given by 
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7 = \J e 2 nj e^mc 2 . (15) 

For copper e.g. 7" 1 ~ 1.8 x 10 -8 m (in the theory of superconductivity 7 _1 is known as 
penetration depth 8 ). The last expression in eq.(14), corresponding to the case in which 
7 2 ai? ^> 1, is independent of geometrical details. The following relation 

(j A ) ~ -(e 2 n/27ri?m)$ B (16) 

has been used in the derivation of eq.(14). Eq.(14) gives an expression for the flux $g in 
terms of p, i.e. in terms of the total angular momentum M = pN. 

For this model the magnetic energy £p , defined in eq.(12), is negative: 

h 2 u 2 2 fi 2 /; 2 Nh 2 n 2 

£ F ~ -(l/2)0 P )$ fl aL = -N^-(l + ^i-)- 1 * -N±±s * (17) 

2mR z j z aR 2mR i 2m 

Second model. - Now let us consider a second model, consisting of a long, thin cylindrical 
conductor, with a current flowing along it. Let L be the length and p the radius, with p Q <C L. 
We assume an uniform stationary canonical current, flowing inside the conductor, given by 

(jp) = —em~ 1 nhp, (18) 

where hp = N~ 1 hP is the canonical momentum per electron (in this case the total number 
of electrons is given by N — irp^Ln). 

For a very long conductor we may assume that A c is approximately oriented in the direction 
of the current (z-direction) and that it depends, approximately, on the distance p from the 
axis of the cylinder only. Putting x = jp, we obtain from eq.(6) 

V 2 A C - 6{x Q ~ x)A c = 6{x - x)he~ l p. (19) 
The solution to eq.(19) has the form 

A c = ~he~ x p + 9(x - x)A mt + 6{x - x )A exU (20) 
where Ai nt is solution to 



Ut +-^-A int =0 (21) 



dx 2 x dx 



and A ext is solution to 



d 2 A ext 1 dA ext _ ^ 
dx 2 x dx 

We assume that, for p -C p -C L, the solution of eq.(22) must coincide with the expression 
for the potential generated by a current ix flowing through a one-dimensional wire 



mp) * 2 ^2 fy^p 2 + z2 y 1/2 ( 23 ) 

* l n |iL + ( l + ^)i/ 2 |^_^_l n ^. 



2-KtQC 2 2p Ap 2 27re c 2 

Taking into account the regularity of A c for x — 0, as well as the continuity of both A c and 
its gradient across the surface x — xq, we obtain 



F. Miglietta : Magnetic energy of a quantum current 



5 




0.5 1 1.5 2 2.5 3 u 

Fig. 1 - Plot of eF versus xq, for a conducting wire .The function eF represents the ratio between the 
total current %t and the canonical one i p .The reduced radius xq = 7P0 represents the radius of the 
wire in units with 7 -1 = 1, where 7 -1 is the penetration depth. For copper e.g. xq = 1 corresponds 
to a radius po — 1.8 x 10 -8 m. Corrections due to surface effects can be expected. The plot is given 
for L/po = 10 4 , where L is the length and po the radius of the conducting wire. 



A int = -[iT/2ne c 2 xoIi(x )] Iq(x), (24) 

where the I's are modified Bessel functions (imaginary-argument Bessel functions) We 
obtain also 



27re c 2 x n x I 1 (x ) 
By a comparison of eq.(20), for x > xq, with eq.(23) we obtain 



A ^ = -^ZT^l ln — + I t TZ J - ( 25 ) 



i T = e F (x )i p , (26) 

with i p = TrpQ(jp) and ep given by 

e F = 2[x 2 ln(X L /x )+x () Io(xo)/Ii{x )}- 1 , (27) 

where we have put Xr, — -fL. Physically eq.(26) is a consequence of the skin-effect. A plot 
for ep versus xq is given in Fig.l, for Xl/xq = 10 4 . 

According to eq.(12) the magnetic energy £p is given by 

E F ~ -(Nh 2 p 2 /2m)[l - e F {x )} ^ -(Nh 2 p 2 /2m). (28) 

The last expression in eq.(28) represents the asymptotic limit for xq 1 (bulk limit). It 
is independent of geometric details and coincides with the r.h.s. of eq.(17) obtained for the 
previous model. Such a coincidence ascribes a sort of universality to this result. In the 
derivation of eq.(28) the following identity 

dx xlo(x) — xqIi (xq) (29) 



has been used. 
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Conclusions. Let £m be the ground-state energy for a system of electrons, in the 
absence of current. Let A£m be the variation for £m, corresponding to a translation in 
momentum space, assuming that such a translation is allowed by the electron state. According 
to eq.(12), the corresponding variation for the total energy of the system is given by 



A£ = £ F + A£ M . (30) 

Since the magnetic energy £f is negative, a physical situation, in which A£ is negative, is 
possible. In such a case the true ground state of the system would be characterized by a 
spontaneous current, as supposed in Ref.s |1| and [Sl- 
it can be demonstrated that, as a consequence of eq.(28), the ground state of a system 
of electrons in a periodic potential, with periodic boundary conditions, is at least degenerate, 
up to finite-size effects. In fact, let \i/jo) be the ground state of Hm and let £m be the 
corresponding eigenvalue. In the state li/'o) the total canonical momentum, as well as the 
total current, vanish. Let l^') be the state obtained from \ipa) through a rigid translation in 
momentum space 

W)=exp\in- 1 p-J2r j ] |Vo>. (31) 

We obtain 



W\HmW) = £m + N{h 2 p 2 /2m), (32) 

and, according to eq.(28), A£ ~ 0. However this result, since it holds up to finite-size effects 
only, is not sufficient to allow a spontaneous current. 

The possibility for A£ < is investigated, for the sake of simplicity, in the case of a 
system of Bloch electrons in a conducting band. However a similar analysis could be applied, 
in principle, to more complicated electron states. In this simple case we have 

£m = {V/8tt 3 ) J alkn{k)e(k) (33) 

and [HI 



p = {V/8ir 3 Nh) J dkn{k){4> k \{-ihV)\^ k ) = (mV/8ir 3 Nh 2 ) J dkn{k)V k e{k). (34) 

Neglecting the magnetic energy of eq.(28), in the ground state one would obtain, for n(k), the 
zero-temperature Fermi-Dirac distribution function 



n F (k) = n F (e(k)) = 26{e F - e(fc)) (35) 

and correspondingly p = and vanishing current. However, owing to the negative value of £ F , 
it may happen that the solution given in eq.(35) is unstable, when the magnetic interaction 
is switched on. In fact let us consider the following distribution function 



n(k) = n F (k — e z q), (36) 

obtained from n F through a shift of the argument. We observe that the transformation of 
the distribution function given in eq.(36) represents a sort of energy-weighted translation in 
momentum space, not a rigid one as given in eq.(31). In this case we have 
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P = I dk ^(^) 2 = -^p- 2 I dk 5(e - e F ){^) 2 (37) 

S^NTi 2 J de y dk z ' Air 3 Nh 2 J { ' dk/ 

and 

AP Vq 2 f dn F de a Nh 2 

where the parity of e(k) has been used. Now eq.(30) reads 

AS ~ (Nh 2 /2m)p(q- p) 
and the condition for A£ < consists of p > q, i. e. 

(V/^N) J dkS(e-e F ){de/dk z ) 2 > (h 2 /m), 

where m is the physical mass of the electron. 

In the simple case of a spherically symmetric e we obtain from eq.(40) 

kp l \de/dk\ k=kF > (h 2 /m) : (41) 

where the expression k F = S^NV' 1 has been used for the Fermi momentum kp. This is 
the condition for a superconductivity a la Heisenberg 4 . Finite-size effects could destroy the 
superconductivity e.g. in the case of a very small ring, according to eq.(17), or of a very thin 
wire, according to eq.(28). 

For a free-electron gas in the bulk limit, the l.h.s. and the r.h.s. of eq.(41) coincide. This 
result indicates that, for such a system, the ground state is degenerate, up to finite-size effects, 
provided that the magnetic self-interaction energy is taken into account properly. 



(38) 



(39) 
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